Gravitino production in the warm inflationary scenario 
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We estimate the production of gravitinos during and after the end of a period of warm inflation, a 
model in which radiation is produced continuously as the field rolls down the potential producing 
dissipation. We find that gravitino production is efficient for models in the strong dissipation 
regime, with the result that standard nucleosynthesis is disrupted unless the magnitude of the 
inflaton potential is very small. Combining this with the constraint from the thermal production of 
adiabatic density perturbations we find the dissipation rate must be extraordinarily strong, or that 
the potential is very flat. 
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I. INTRODUCTION 

The warm inflation scenario Q| is an unusual vari- 
ant on the inflationary cosmology, in which the inflaton 
has significant interactions during the inflationary epoch 
leading to a continuous production of radiation. The 
backreaction of this production on the inflaton field ap- 
pears as a viscosity, slowing down the scalar field evo- 
lution and hence aiding slow-roll inflation 0. In such 
a scenario, inflation can proceed with potentials steeper 
than those in standard chaotic scenarios. 

Issues concerning the implementation of warm infla- 
tion within a realistic particle physics context have yet 
to be studied to the same depth as the standard infla- 
tionary scenario HQ . Nevertheless, given that the warm 
inflationary scenario is very different phenomenologically 
from the usual picture, it makes good sense to examine 
the extent to which its phenomenology is consistent with 
observations. The two main purposes of inflation are to 
provide a large, nearly homogeneous, patch in the Uni- 
verse within which structure formation can take place, 
and to ensure that unwanted relic particles do not spoil 
the successes of the standard hot big bang cosmology. 
The first of these has seen a reasonable amount of study 
||, and so we will consider an example of the latter. 

In the context of modern particle physics, the most 
troublesome relics are the gravitino and the moduli fields 
||. We will consider the gravitino, whose existence arises 
as the supersymmetric partner of the graviton, and whose 
mass is expected to be order of 1 TeV. It is a cosmo- 
logical threat because if produced in enough abundance 
in the early Universe, it is sufficiently long lived to sur- 
vive until after nucleosynthesis, at which point its decays 
spoil the element abundances . To avoid this, the ratio 
of gravitino to photon number densities must be below 
about 10~ 12 . The gravitino may be produced both by 
interactions within a thermal bath and by various 
non-thermal processes ||. In conventional inflationary 
scenarios, the former gives an important upper limit on 
the reheat temperature, while the latter may constrain 
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many possible physical processes. 

In this paper we explore the consequences of gravitino 
production during and after warm inflation. Warm infla- 
tion differs from conventional inflation in that radiation is 
constantly produced during inflation, and the radiation 
density decreases monotonically throughout the evolu- 
tion, with inflation ending when the radiation density 
overtakes the inflaton energy density. There is therefore 
continuous gravitino production during inflation, and 
also no delay in post-inflationary thermal production due 
to an intervening (p)reheating period. Consequently the 
gravitino bound is much harder to satisfy. We will show 
that the abundance of gravitinos produced during infla- 
tion is similar to that produced after inflation, and assess 
the strength of the constraints this imposes on warm in- 
flation model-building. 



II. EVOLUTION OF FIELDS DURING WARM 
INFLATION 

A. Dissipation effects during inflation 

We review the dynamics of warm inflation closely fol- 
lowing Taylor and Berera ||, where full details can be 
found. Warm inflation is distinguished from ordinary in- 
flation by the presence of a viscous damping during the 
inflationary evolution, so that the inflaton field <p satisfies 
the equation 



4> + (3H + r)0 + V = . 



(1) 



Here H = a/a is the Hubble parameter, a is the cosmo- 
logical expansion factor, and T is the dissipation coeffi- 
cient. V((f>) is the potential of the inflationary field. For 
simplicity we assume a spatially-flat universe through- 
out. A dot denotes differentiation with respect to time 
and a prime with respect to <fi. 

The energy density of relativistic species, p ra d, follows 
from energy conservation as 



1 



p rad + 4_ffp rad = T(f> 



(2) 



These equations are completed by the Fricdmann equa- 
tion 



H 2 = 



8tt 
3m P1 



{p<t> + Prad) 



(3) 



where mpi is the Planck mass and the energy density of 
the inflaton field is 

P* = V{<t>) + \tf. (4) 

During an inflationary era the potential field dominates 
both the kinetic energy of the inflationary field and the 
energy density of the radiation, so the Friedmann equa- 
tion can be reduced to 

Assuming the slow-roll condition, the equa- 

tion of motion for the inflaton reduces to 

: V 



where 



3JT(l + r) 

r 

= 3H ' 



(6) 



(7) 



is a dimensionless dissipation coefficient, whose value in 
general is time dependent. During the inflationary period 
the production of radiation will soon settle into a stable 
state, where p ra d <C r</> 2 , giving 



Prad 



(8) 
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where g* is the number of relativistic degrees of freedom. 
In the Standard Model of particle physics = 106.75 
for T > 300 GeV, while in the minimal supersymmetric 
model (MSSM) we shall be assuming here this rises to 
<7* = 228.75 once the temperature is above the mass of 
the supersymmetric particles. 

Combining Eqs. ([§, (§) and (§) wc find 
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where 
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(10) 



is the usual inflationary slow-roll parameter. Eq. (|^) 
holds for all values of r. In the limit r — > dissipation is 
switched off, and the radiation field vanishes as 



Prad — — / 



(11) 



In this regime, increasing the dissipation factor r in- 
creases the decay of the inflaton field into radiation, while 
having no effect on the evolution of the inflaton field. For 
fixed e the fractional density of radiation is at a maximum 
when r = 1 (equivalently T = 3H). 



B. The strong dissipation regime 

Our main focus will be the strong dissipation regime, 
where the differences from standard inflation are most 
pronounced. In the regime of strong dissipation r 1 
the radiation field is given by || 



Prad 
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(12) 



The major effect of increasing the dissipation factor is to 
heavily dampen the evolution of the inflaton field, slowing 
its evolution down the potential and decreasing the decay 
into radiation. 

The conditions for slow-roll and warm inflation to oc- 
cur are 

e<2r, (13) 
and for an extended period of inflation we need 



where 
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Eqs. ( |l3| ) and (|1J) relax the usual constraints on the infla- 
tionary potential. Supercooled inflation ends when e«l, 
while warm inflation takes place until 



2r, 



(16) 



when the radiation energy density starts to dominate the 
energy density of the inflaton field. At this point the 
universe makes a smooth transition from the inflation 
phase to a radiation-dominated, hot Friedmann model. 

To illustrate the evolution of warm inflation we will 
consider polynomial potentials of the form 



V{4>) = Am 4 X. 

m 



(17) 



where we allow a to be a positive real number. It is 
important to note that the only fully complete model of 
warm inflation, where a period of warm inflation comes 
to a natural end with the field finishing in the minimum 
of its potential, occurs when a = 2. For other positive 
even integers the potential has a suitable minimum at 
the origin. However for a > 4 the radiation density falls 
compared to the inflaton energy and so warm inflation 
does not take hold, as can be seen from Eq. ( |l2| ) since 



and 



The case a = 4 is spe- 



cial in that the densities of the two components remain 
in fixed proportion; inflation proceeds forever with the 
dissipation easing the inflaton asymptotically into the 
minimum. Modification to the potential would be re- 
quired to end inflation. We will also consider odd and 
non-integral values of a (taking the field to have positive 
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III. THERMAL PRODUCTION OF GRAVITINOS 




FIG. 1. The evolution of energy density of the inflaton and 
radiation fields as a function of number of e-folds from the end 
of warm inflation, for the inflationary potential V = m 2 (f> 2 /2 
(see text for details). The energy density of the inflationary 
field (upper curves) and the radiation field (lower curves) are 
taken as initially equal, and the vertical scale is arbitrary. 



values). Such potentials have no minimum, and indeed 
may be ill-defined for negative <j>; we include them merely 
to illustrate the effects of modifications to the slope of the 
inflaton potential. Our range of investigation will cover 
1 < a < 4. 

For the polynomial potential the slow-roll parameter 
e = cv 2 mpj/167T0 2 , and the number of e-folds of expansion 
to the end of inflation, N(<j>), is § 



N((j>) = In 



^cnd 



a p r ad 



(18) 



where and p ra( j are the energy densities near the start 
of warm inflation once stable radiation production has 
been established. If warm inflation starts with an initial 
stable ratio p^/p T3 A, it will take N ~ p^/p Ia d e-folds 
before the radiation and the vacuum energy are equal. 

Figure 1 shows the evolution of the energy densities 
of the inflaton and radiation fields during an inflationary 
era, when a = 2. We numerically solved Eqs. (Q), (0) and 
(§) as in Ref. 0. The number of e-folds is N = 100, and 
we began the model with p ra( j = p<j>, although the evolu- 
tion is insensitive to the initial conditions and soon settles 
into its stable configuration. The choice of parameters 
was made so that the amplitude of thermally-produced 
adiabatic perturbations generated during warm inflation 
agrees with the amplitude of temperature fluctuations 
in the CMB measured by COBE, fe = 2x 10~ 5 (see 
Section IV.B and Ref. [3 for details), with T 

4 m P1 . 



m = 10- 8 m P i, and V" 1 / 4 - 10~ 4 - 
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Gravitinos are too weakly interacting to be able to 
reach thermal equilibrium with a radiation bath unless 
the temperature is around the Planck temperature. How- 
ever, although interactions are negligible once they form, 
they can be created by two-body processes such as 



7 + 7 -> 9 + 9 ■ 



(19) 



The single-particle decay rate for gravitinos gives a life- 



time of order rn^ 2 



/m 2 



pi- 



A. Thermal production during warm inflation 

The number density, TI3/2, of gravitinos produced from 
the thermal bath is governed by the equation || 



«3/a +3-H%/ 2 = (o- 3 / 2 |w|)n 



rad ! 



where 



™rad 



C(3) 



g*T 3 



n 98« 1 / 4 n 3 / 4 
U - 28 3* Prad' 



(20) 



(21) 



is the number density of particles in the thermal bath, 
and C(3) = 1.202. Provided the typical particle energies 
well exceed the gravitino mass, which is always an excel- 
lent approximation for us, the effective total cross-section 
in a thermal bath for gravitino production, including all 



particle channels, has been computed as a^ 2 



250/toIj 



|)8|,p.0|. For convenience we have defined a mean creation 
rate per particle species as 03/2 ~ 250/ i g 2 rn| 1 , the square 
arising as the production is a two-body process]^] We will 
be considering temperatures hot enough that all super- 
symmetric species participate in the thermal bath. The 
factor v « 1 is the velocity of the produced gravitinos. 
The gravitino production is small enough that backre- 
action on the radiation density can be neglected, as can 
reactions destroying gravitinos by interactions or decays. 

After a short period of inflation the gravitino produc 
tion rate becomes stable, 71.3/2 -C (c3/2M)^r 
gravitino number density is given by 



2 ad , and the 



n 3/2 — fg^rad ; 



where 



(03/2M) 
3H 



rirad 1 



(22) 



(23) 



is the dimcnsionless production rate in units of the Hub- 
ble expansion. We define the yield of gravitinos as 



*Our notation differs from that of Ref. [g], who define n ra( j 
to be the number density of a single degree of freedom and 
use the total cross-section. 
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FIG. 2. The evolution of number density of the radiation 
and gravitinos as a function of number of e-folds from the end 
of warm inflation, for the inflationary potential V = m 2 (f> 2 /2. 
Number densities are expressed in units of the radiation num- 
ber density at the 50 th e-fold. 



FIG. 3. The evolution of gravitino yield for a = 2, 3.8 and 
4. The thick lines are calculated numerically, while the lighter 
lines are from the analytic expression Eq. (^). The warm 
inflationary parameters are chosen to produce the observed 
amplitude of adiabatic density perturbations (see text). 



) 



3/2 



n 3 /2 



(24) 



where the entropy s is equal to 3.6n ra( j in the high-energy 
regime. Hence during warm inflation the yield is simply 



Y 3/2 = 



3.6 



(25) 



As long as the stable production hypothesis is valid, the 
ultimate yield will only depend on the situation at the 
end of warm inflation. 

Figure ^ shows the evolution of the number densities 
of the radiation and gravitino populations during warm 
inflation. The number densities are normalised to the 
number density of radiation at the 50 th e-fold, and are 
typically ~ 10 12 to p1 , and n 3 / 2 ~ 10 20 mp,. The 
choice of model parameters is the same as for Figure 1. 

The gravitino yield is given by Eqs. @, (|||), ((25|) , 
© and (|) 



(26) 



For a fixed potential, the yield during warm inflation in- 
creases with increasing potential magnitude and slope, as 
there is greater dissipation into radiation, while the yield 
decreases with increasing dissipation factor as strong dis- 
sipation will dampen the decay process. 

For polynomial potentials of the form of Eq. ([l7]) the 
yield is given by 



3/2 



mpi v r / 



m 



(5a-12)/g 



(27) 



For this type of potential we see that a = 12/5 is a critical 
slope for gravitino production, leading to a constant yield 
as a function of time. Expressing Eq. ( p7| ) in terms of 
the number of e-folds until the end of warm inflation 



Y, 



3/2 



(4 - a) 



(5a-12)/4(4-a) 



(28) 



For a < 12/5 the yield increases as a function of time, 
or e-folds, with a maximum at the end of warm inflation. 
For a > 12/5 the yield is a decreasing function of time. 

Figure g shows the evolution of the yield Y 3 / 2 dur- 
ing the warm inflation phase, for polynomial potentials 
with a = 2, 3.8 and 4. The solid lines are calculated 
numerically, while the lighter lines are from the analytic 
expression, Eq. (|27]). Since the gravitino number den- 
sity quickly settles into its stable production state, the 
analytic expression accurately describes the evolution. 

As the slope is increased to a — ► 4, the timescale for 
inflatio n is a symptotically stretched out, as discussed in 
Section |ll B[ The dependence on e-folds asymptotically 



becomes Y 3 / 2 ~ e^ 4 = (a/a e nd) -1 • This weak depen- 
dence can be seen in Figure 3. 

At the end of warm inflation e = 2r and the yield is 
given by 



4 



Y m-^ Vi — (29) 

Hence at the end of warm inflation the yield only de- 
pends on the magnitude of the potential. The end yield 
is directly related to the final temperature, T en( j, by 

F 3/2 (T cnd )^1.5 5 7 3/2 (^) . (30) 



B. Evolution of gravitinos after warm inflation 



constrain the warm inflationary parameters. In standard 
inflation, there are essentially two free parameters, the 
amplitude of the inflaton potential, V, and its slope, e. 
In general these can be constrained by the gravitino yield 
and the amplitude of adiabatic perturbations, the lat- 
ter being constrained by the observed fluctuations in the 
microwave background. In the following section we shall 
show that the main constraint from the gravitino produc- 
tion is on the magnitude of the inflaton potential, while 
the adiabatic density perturbations constrain the dissi- 
pation factor r. We begin with the constraint from the 
gravitino production. 



Gravitinos have a sufficiently long decay time that 
they survive beyond nucleosynthesis. Their decays de- 
stroy 4 He and D nuclei by photodissociation, and if the 
gravitino abundance is high enough this will disrupt the 
successful predictions of nucleosynthesis. The important 
quantity is the ratio of the gravitino to entropy densities. 

The entropy is a particularly useful quantity to follow, 
as the comoving entropy is conserved not only during 
normal expansion but also during epochs where species 
fall out of thermal equilibrium and annihilate, changing 
the number of particle species in the thermal bath, so 
s oc I/a 3 always.^ The evolution of gravitinos, Eq. (20), 
can be rewritten 



"■3/2 + 3^713/2 = sY 3/2 = (cr 3 /2 M>™?ad 



(31) 



The yield produced once radiation domination begins is 
readily calculated using the normal radiation-dominated 
solution. It is dominated by early time production, as 
is well known, and the total yield from the beginning of 
radiation domination is 



^3/2 - Y 



end 

3/2 



39 
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m P i 



i 1/3 



(Tend) 



i 7 d 3 ( 



(T) . (32) 



Once T -C T cn d, the right-hand side of Eq. (|32j) equals 
the yield at the end of warm inflation, given by Eq. (29), 
since at the end of warm inflation p$ ~ p ra d- We con- 
clude therefore that the production of gravitinos during 
warm inflation can only lead to at most a factor of two 
enhancement over the total production at the end of the 
warm inflation era. 



IV. CONSTRAINTS ON WARM INFLATION 

Having calculated the yield of gravitinos during the 
warm inflationary era, and shown that the yield does 
not significantly change afterwards, we now use this to 



^An exception to this would be if there were particles with 
late out-of-equilibrium decays. 



A. The constraint from nucleosynthesis 

Avoiding overproduction of D + 3 He constrains the 
ratio of gravitinos to photons at the end of warm infla- 
tion. The details of the constraint depend on the grav- 
itino mass but for our purposes we can safely adopt 
a conservative limit 



^3/2 < 10 



-12 



(33) 



Combining this with Eq. (B9[ ) for the yield at the end of 
warm inflation, and Eq. (|32^ , for the subsequent produc- 
tion during radiation domination, we find the following 
constraint on the magnitude of the inflation potential 



V 1 ' 4 <2xHT 13 gl /4 m Ph 



(34) 



This can be expressed in terms of the temperature at the 
end of the warm inflationary phase, via Eq. ( |30| ) , giving 



T cn d < 8 x 10 6 5 * /2 GeV . 



(35) 



With = 228.75 for the MSSM, this gives a constraint 
in good agreement with the standard result ||, though 
slightly weaker due to our adoption of a conservative con- 
straint on Y3/2. This is as expected since warm inflation 
has not greatly enhanced the gravitino yield. The main 
difference in warm inflation is how this constraint inter- 
acts with other constraints on the scenario. 



B. The constraint from the amplitude of adiabatic 
density perturbations 

In addition to the yield we can add the independent 
constraint on warm inflation parameters from the ampli- 
tude of perturbations, S^, produced from thermal fluctu- 
ations during the warm inflationary era ^j, 



5 2 h 



.57,- ^ (3 
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'PI, 
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(36) 



The observational constraint from adiabatic density per- 
turbations does not actually constrain the amplitude of 
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the potential field alone [^), since the amplitude is ob- 
served at a fixed e-fold from the end of warm inflation, 
but rather the combination r 2 V = r 2 m P1 /247r. As we 
usually assume the amplitude of perturbations is mea- 
sured at the 50 th e-fold from the end of warm inflation, 
some scaling has to be made to match the constraint from 
the yield at the end of warm inflation. We may write the 
number of e-folds from the end of warm inflation as 



N ~ A 



Gr 1 )' 



(37) 



where A depends on the shape of the potential. For poly- 
nomial potentials 



.4 



2a 



(4 -a) 



(38) 



The amplitude of perturbations then constrains the dis- 
sipation factor 



10" 



4 x 10" 
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(39) 



As r = T/H and H ~ V 1 / 2 , we can combine Eq. j39| ) 
with the constraint on the gravitino yield, giving a con- 
straint on the dimensionless dissipation rate 



r > 2 x 10 1 



4 x 10" 10 



2/3 



-1/2 



(40) 



This equation is our main result. It shows that if warm 
inflation is to simultaneously produce density perturba- 
tions of a satisfactory magnitude and avoid overgenerat- 
ing gravitinos, then the dimensionless dissipation must 
be extremely high. The main reason for this is that high 
dissipation increases the magnitude of density perturba- 
tions, allowing the inflationary energy scale to be normal- 
ized down. With sufficient dissipation, the energy scale 
becomes low enough that the gravitino yield is sufficiently 
suppressed. 

One caveat to Eq. ([hj|) is the factor A, relating the 
ratio r/e to the number of e-folds before the end of warm 
inflation, when the density perturbations where formed. 
For polynomial potentials this is typically of order unity, 
but for potentials with very flat slopes, a — > 0, the dis- 
sipation rate can be arbitrarily small, as the radiation 
production is suppressed. 



V. CONCLUSIONS 

A crucial role of inflation is to ensure that unwanted 
relic particles do not survive with abundances capable 
of spoiling the standard hot big bang model. We have 
studied the production of gravitinos during and after a 



warm inflationary era, and combined the constraint this 
gives with the requirement that density perturbations of 
the correct magnitude are generated. 

We have found that although there is continuous grav- 
itino production during warm inflation from interactions 
in the thermal bath, this does not in itself lead to a 
very significant extra yield of gravitinos over and above 
that produced at the end of warm inflation. Neverthe- 
less, avoiding overproduction of gravitinos is much more 
challenging than in conventional inflationary scenarios, 
because the radiation density is monotonically decreas- 
ing throughout the evolution. Satisfying the gravitino 
bound requires that the potential energy at the end of 
warm inflation be very small, U 1 / 4 < 10~ 9 mpi. In or- 
der for density perturbations to have the correct mag- 
nitude requires a dissipation factor T w 10~ 6 mpi, and 
hence a dimensionless dissipation r > 10 12 unless the 
slope of the potential is extremely flat. Either way, it is 
clear that evading overproduction of gravitinos strongly 
constrains the warm inflation scenario, requiring dimen- 
sionless numbers many orders of magnitude away from 
unity. 
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